Spin-dependent dipole excitation in alkali-metal nanoparticles 
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We study the spin-dependent electronic excitations in alkali-metal nanoparticles. Using numerical 
and analytical approaches, we focus on the resonances in the response to spin-dependent dipole 
fields. In the spin-dipole absorption spectrum for closed-shell systems, we investigate in detail the 
lowest-energy excitation, the "surface paramagnon" predicted by L. Serra et al. [Phys. Rev. A 47, 
R1601 (1993)]. We estimate its frequency from simple assumptions for the dynamical magnetization 
density. In addition, we numerically determine the dynamical magnetization density for all low- 
energy spin-dipole modes in the spectrum. Those many-body excitations can be traced back to 
particle-hole excitations of the noninteracting system. Thus, we argue that the spin-dipole modes 
are not of collective nature. In open-shell systems, the spin-dipole response to an electrical dipole 
field is found to increase proportionally with the ground-state spin polarization. 

PACS numbers: 78.67.Bf, 73.21.-b, 73.22.-f 
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I. INTRODUCTION 

The optical absorption of small metal particles is 
dominated by the surface plasmon resonance^ 2 -^^. In 
this collective excitation, the center-of-mass of the elec- 
trons moves back and forth with respect to the posi- 
tive background leading to an oscillating charge dipole. 
Pump-probe experiments using femtosecond laser pulses 
have been widely used to study the relaxation of these 
excitations^. Within the surface plasmon excitation all 
the electrons oscillate in phase, irrespective of their spin. 
Thus, its study only yields information on the charge dy- 
namics. In order to address the spin dynamics of nano- 
objects, time-resolved magneto-optical Kerr effect mea- 
surements have been performed recently, yielding the full 
trajectory of the magnetization in real space for optically 
excited superparamagnetic nanoparticles^. 

A substantial interest in the theoretical description of 
spin dynamics in nano-objects was aroused by the work 
of Serra et who found that a strong peak in the spin 
dipole absorption spectrum of alkaline nanoparticles ex- 
hausts a large fraction of the energy-weighted sum rule 
and drew the conclusion that this peak corresponds to a 
collective spin mode. In this excitation named "surface 
paramagnon" , the spin degree of freedom appears in a 
crucial manner. While the surface plasmon can be ex- 
cited by a dipole electric field, the surface paramagnon 
results when an excitation acts differently for spin up 
and spin down electrons. Such a field can be realized 
experimentally through the magnetic field component of 
electromagnetic waves having a wavelength considerably 
longer than the size of the system^. 

For the case of spherically symmetric systems with 
zero total spin, in the spin-dipole mode the center-of- 
mass of the electron system does not move. Therefore, 
such an excitation does not couple to an electric field in 
closed-shell systems. Since transitions induced by elec- 
tric dipole fields are the dominant mechanism, the sur- 



face paramagnon is difficult to observe in the optical ab- 
sorption spectrum of such nanoparticles. However, for 
nanoparticles with no spherical symmetry and/or open- 
shell electronic systems, the spin dipole couples to elec- 
tric dipole field s 10 ' 11 ! 12 . For relatively small systems, 
it has been concluded from calculations of the time- 
evolution of strong excitations that the coupling between 
the charge and the spin-modes is not crucially modified 
when the excitation strength is increased into the non- 
linear regime^ -. In systems that lack spherical symmetry 
scissor modes exist which can also be coupled to the spin- 
modes^. 

The difficulty to detect the spin-dipole mode in ro- 
tationally invariant nano-objects becomes less restric- 
tive when one studies semiconductor quantum dots in- 
stead of metallic nanoparticles. For the typical sizes and 
electronic densities of quantum dots the energies of the 
charge and spin dipole excitations are of the order of 
a few meV. Therefore Raman scattering of visible light 
can be used (with selection rules depending on the po- 
larization geometryii) to detect and study both kinds of 
excitations^. Finally, the strong electronic confinement 
of disc-shaped quantum dots in semiconductor hetero- 
structures results in very sharp resonances for dipolar 
excitations. The enhanced damping of the lowest spin- 
dipole mode in the presence of a weak magnetic field that 
splits the single-particle excitations has been taken as an 
indication of the collective character of the lowest spin- 
dipole modei^. 

The availability of experimental data has motivated 
theoretical work on the charge and spin density excita- 
tions in semiconductor few-electron quantum dots (see 
for example Ref. Il6l). In particular, the induced magne- 
tization density has been studied^ and a correspondence 
between the spin-dipole modes and single-particle excita- 
tions has been observed numerically^. The small devi- 
ation of the spin density resonance energies from single- 
particle excitation energies has been explained by the 
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absence of long-range Coulomb interaction terms in the 
energy of the spin modes. In addition, the excitation en- 
ergies have been studied using the Quantum Monte Carlo 
technique^ and within a semicalssical approach for the 
time-dependent charge and spin density oscillations^. 

In this work, we use analytical (mean-field) and nu- 
merical time-dependent local spin-density approximation 
(TDLSDA) approaches to study in detail the physics of 
the spin-dipole modes in alkali-metal particles. We ad- 
dress the important questions related to the specificity 
of the lowest frequency resonance as compared with the 
other excitations by identifying the modes with the cor- 
responding dynamical magnetization densities. 

By studying the evolution of the spin modes with the 
interaction we provide clear-cut arguments in the discus- 
sion over the collective versus single-particle nature of 
the surface paramagnon. In addition we obtain the size 
scaling of the lowest resonance frequency and relate it 
with spill-out effects. In the case of open-shell systems 
we analyze the nature of the different spin-dipole excita- 
tions and the relationship with the charge modes. 

While our numerical results are worked out for the case 
of alkaline metal nanoparticles, most of our general con- 
clusions hold for a broad class of nanosystems, including 
semiconductor quantum dots. 

The paper is organized as follows. In Sec. UH we 
describe the model for the electron dynamics in metal 
nanoparticles and the numerical method we employ. We 
also present the computed spin-dipole absorption spec- 
trum and the corresponding dynamical magnetization 
density for a typical example of a closed-shell system. 
In Sec. IIII1 we use a phenomenological approach to de- 
scribe the energetically lowest spin dipole excitation and 
derive results for its frequency based on plausible as- 
sumptions on the dynamical magnetization density. We 
compare the numerically obtained frequencies with the 
phenomenological ones resulting from increasingly accu- 
rate descriptions of the electronic dynamics. In Sec. IIV1 
we present numerical results for the full absorption cross 
section in the frequency regime below the surface plas- 
mon frequency. We follow the evolution of the absorp- 
tion spectrum with the strength of electron-electron in- 
teractions and find a one-to-one correspondence of the 
spin-dipole modes with the particle- hole excitations. In 
Sec. El we study open shell clusters and discuss the pos- 
sibility to observe the spin-dipole modes in the electric 
dipole absorption spectrum. We provide our conclusions 
in Sec. IVII In the appendices we present the details of 
our LSDA parametrization and the calculations for the 
case of a non-uniform ground-state electron density. 



II. NUMERICAL APPROACH TO 
SPIN-DIPOLE EXCITATIONS 

In our study of the electronic excitations of nanoparti- 
cles, we restrict ourselves to the electronic degrees of free- 
dom and describe the confining effect of the ionic back- 
ground by a spherical jellium model with sharp bound- 



aries. Such a simplification can be justified for not too 
small metal particles. Furthermore, we do not consider 
thermal effects and therefore choose to work at zero tem- 
perature. 

We start by introducing the formalism underlying 
the numerical approach to the absorption cross sec- 
tion corresponding to spin-dependent excitations of our 
model nanoparticles. We follow the formulation of the 
TDLSDA 2 - as it is presented in Ref. |lf| In this frame- 
work the electronic system is described in atomic units 
(h — to — e = 47reo = 1) using the Kohn-Sham equations 
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where <f>Z is the fe th Kohn-Sham wave-function with the 
quantum number a = {t,4} describing spin projection 
onto the z-axis. These wave-functions allow us to define 
the spin-dependent electron densities 



n a (r,t) 
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where the sum runs over the single-particle like Kohn 
Sham levels k that contribute to the many-body density. 
The electron density, magnetization density, and spin po- 
larization are obtained from n a , respectively, as 



n = rv + , 
to = ?i 1 — n l , 
£ = m/n. 



(3a) 
(3b) 
(3c) 



The effective potential in the Kohn-Sham equations 
can be written as 

V&(r, t) = V c (r) + V u (r, t) + V£(r, t) + V£(r, t) , (4) 

where V c represents the confinement due to the jel- 
lium background, Vh is the Hartree potential, V£ c is the 
exchange-correlation potential and stands for the ex- 
ternal perturbation. The local character of the approxi- 
mation is reflected by the choice 
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where e xc (n^ , n^-) stands for the exchange-correlation en- 
ergy density for which we use the parametrization of 
Perdew and Zunger 21 reproduced in App. [A"l 

Within linear response theory we write the density 
changes induced by the external perturbation as 



(6) 



where V^(r,w) is the Fourier transform of the time- 
dependent external potential, and the response functions 
X <T<T obey the Dyson equation 
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X aa '(r,r',u;) = X ™'(r,r',u;) + a c J2 // dridr 2 ^(r.n.w) 



ri - r 2 
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where we introduced the parameter a c = 1 which will later allow us to modulate artificially the importance of the 
electron-electron interactions in model calculations. The kernel of Eq. [7] is given by 
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The non-interacting response function is diagonal in the spin indices and given by the density-density correlator 

1 1 1 
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which can be expressed in terms of the retarded Green functions^^. We have chosen the imaginary part in the 
denominator that ensures the convergence as rj = 8meV (= 2.94 x 1CP 4 in atomic units). The robustness of the final 
result with respect to variations of this parameter has been checked. 

Defining the spin-independent part of the external perturbation V cx . n — (Vj x + V cx )/2 and its spin-dependent 
counterpart V CXtm = (VJ X — Vj x )/2 > the response of the charge and magnetization densities n and m can be expressed 
in matrix form as 



5n(r,w) \ [ dr > f Xnn(r,r',w) Xnm(r,r',w) \ / K x , n (r',w) 
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where the cross-correlations of the charge and spin chan- 
nels are given by 



Xnn/nm = X^ ± X ' 1 + X 1 ' ± X 

X - - -\/TT -|— T/Ti- — -\/J-T zn. v"^^ 
mn/mm A A A \ A 
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Electromagnetic radiation with wavelength much 
larger than the size of the nanoparticles induces dipo- 
lar perturbations. Considering monochromatic light with 
wave- vector k = fcx, linear polarization along y, and 
therefore magnetic field along the z-direction, the dipole 
excitation potentials for charge and spin can be written 
as 



Kx,n(r,i) = F D y, 



(12a) 
(12b) 



with the excitation strengths F n = — £ ma x sin(wi) and 
F m = g[iBE max k cos(wi), the Bohr magneton /ib, and the 
gyromagnetic factor g (3/13 = 1 in atomic units). These 
dipolar perturbations lead to dipolar charge and spin- 
density excitations, and the corresponding polarizabili- 
ties are given by 
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drdr' rr 1 cos 9 cos8' Xab( r , r ', w ) 7 (13) 



with ab={nn, nm, mn, mm}. For the simplicity of nota- 
tions, the polar coordinates (r, 9, ip) are here and hence- 
forth chosen to have the z-axis along the variation of 
the excitation field. This conventional choice allows us 



to treat the charge and spin excitations within the same 
description, but it is not consistent with the example of 
electromagnetic radiation presented above. 

The experimentally relevant quantities are the dipole 
absorption cross-sections 



Sab(w) = — Im [a ab (w) 



(14) 



For spherically symmetric nanoparticles a and S are di- 
agonal in the channel indices a and b. In this case we will 
work with S n = S nn and S m = S mm . The spin-dipole 
absorption spectrum S m (ui) for the closed-shell system 
Na34^2 is shown in the left panel of Fig. [TJ Four peaks 
are observed in the low-energy range below 0.6 mm, with 
the Mie frequency ujm = 3.4 eV which is the classical fre- 
quency of the surface plasmon excitation^. The peak at 
the lowest frequency, labeled (1) in the figure, displays 
the strongest absorption cross-section and corresponds 
to the surface paramagnon described by Serra et aZ£. 
In section iriTl we derive analytical expressions that accu- 
rately describe its frequency and its dependence on the 
size of the nanoparticle. 

The radial part of the magnetization density at reso- 
nance is shown in Fig. [1] (right) for the four peaks ap- 
pearing in the absorption spectrum. The magnetization 
profile for the lowest frequency peak (1) clearly differs 
from the profiles corresponding to the higher frequency 
peaks: it involves considerably stronger magnetization 
densities than the other peaks and, most importantly, 
it displays no significant nodes (except at the center of 
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FIG. 1: Left: spin-dipole absorption spectrum for a Na34 nanoparticle. The frequencies of the horizontal axis are normalized 
to the Mie frequency ujm = 3.4eV (= 0.125 in atomic units). Right: radial part of the magnetization density for each of the 
resonances identified in the left panel, as a function of the radial coordinate (scaled with the radius a of the particle). The 
magnetization profile is scaled with the excitation strength, the mean density and the radial coordinate 



the nanoparticle), whereas the other peaks are associ- 
ated to magnetization profiles with richer node struc- 
tures inside the nanoparticle. For larger systems, even 
more peaks appear, and the corresponding magnetiza- 
tion profiles show more complicated structures with sev- 
eral nodes. However, the particularly simple structure 
of the lowest-frequency peak, and its stronger amplitude, 
persist and thereby point to its special character. 



A. Uniform ground-state density 



The ground-state equilibrium electron density in a 
spherical jellium model of radius a with sharp bound- 
aries can be approximated by a uniform distribution 



III. PHENOMENOLOGICAL APPROACH TO 
SPIN-DIPOLE EXCITATIONS 

The physics of the spin-dipole excitations obtained in 
the previous section can be understood through phe- 
nomenological models. In particular, we present an esti- 
mation of the lowest resonant frequency and compare it 
with results from TDLSDA calculations. 

We consider a spherically symmetric nanoparticle with 
ground-state equilibrium densities nj(r) = np(r) and a 
perturbation in the spin channel such that An^(r) — 
— An'(r). The displacements of the center-of-mass of 
the two spin populations along the z-direction are given 
by 

Z T = -Z l = jLJdrz An T (r), (15) 

where = N/2 is the number of spin-up electrons and 
N the total number of electrons. Since there is no net 
charge displacement, the Hartree term En of Eq. IA5I rc- 
mains unchanged under the perturbation. The changes 
in the other contributions to the total energy can be cal- 
culated from ArJ(r). 



(r) = n^(r) = — 8(a — r) (16) 



inside the sphere, where n = 3N/4ira 3 and 8 denotes the 
Heaviside function. Assuming that the perturbation is a 
dipolar field, the simplest approximation to describe the 
low-energy spin excitations of the system is to postulate 
the tilts 



An T (r) = -An l {r) = -- Q(a-r) (17) 
C 2 



of the spin densities. The characteristic length ( de- 
scribes the magnitude of the excitation. Working in lin- 
ear response, we restrict ourselves to weak excitations 
with ( > o and consider the change in total energy in- 
duced by the above density excitations using the energy 
functionals described in App. [XJ 

The simple form assumed for the spin densities allows 
us to neglect the change of Ek.g ■ To the lowest order in 
the perturbation the other components of the energy are 
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(18a) 



(18b) 



(18c) 



We have defined 



(19) 

with given in App. [£] and r s = (4tt n/3) _1/3 . The 
energy increase due to the spin density displacements 
leads to a restoring force T = —dE/dZ^ and an out- 
of-phase oscillation of the two spin subsystems with a 
frequency 
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We have expressed the result in terms of the classic al 
Mie frequency^ which can be written as ujm = \/N/a 3 — 

(r s ) . 

As compared to the spin dipole, the surface plasmon 
excitation is of quite different nature since it results from 
the oscillation of the total charge. The frequency ujm 
can be obtained following similar lines as those presented 
above, but restricting the restoring force to the Hartree 
contribution. The different nature of the energies in- 
volved in each mode results in a higher frequency for the 
surface plasmon (in the visible part of the spectrum for 
the case of metal nanoparticles) than for the spin dipole 
(in the infrared range). Moreover, ujm is independent of 
the size of the particle, while ojs decreases with the num- 
ber of electrons as A -1 / 3 . T his p ower-law scaling has 
already been obtained in Ref. |23| within sum rule and 
hydrodynamic approaches. It makes the observation of 
the surface paramagnon in large particles more difficult. 

In Fig. we compare the values of u> s d obtained from 
the TDLSDA (filled circles) with the estimate uis of Eq. 
[20] (dotted line). For comparison we also show the numer- 
ically calculated surface plasmon frequencies o; sp (filled 
squares), which approach the classical value ujm for large 
N and display important oscillations for small A [24j |. 
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FIG. 2: Size-dependence of the surface paramagnon fre- 
quency, together with the surface plasmon frequency lj sp 
(filled squares). Filled circles represent ui s d obtained from 
TDLSDA calculations. The dotted line is the estimate ujs 
of Eq. 1201 The open circles stand for the spill-out corrected 
Eq. 1201 where u sp is used instead of uim in the prefactor. The 
pluses depict the semi-analytical result of (|28[) using the mag- 
netization profile arising from Eq. 1261 Inset: Radial variation 
of the ground-state electron density used in Eq. 1261 and ob- 
tained from static LSDA calculations. The dashed, dotted, 
and solid lines are for Na2o, Naio6, and Na306, respectively, n 
is the electron density for bulk Na and the dash-dotted step 
function corresponds to the uniform ionic jellium. 



We can see that the predicted decrease of the spin-dipole 
frequency as iV -1 / 3 is essentially correct. However, Eq. 
[2111 overestimates the actual frequencies. This discrep- 
ancy becomes increasingly important when the size a of 
the nanoparticle diminishes. Two key assumptions in 
the derivation of ([2"0|) become less justified when a gets 
smaller. On one hand, the spill-out effect due to the 
extension of the electron wave-functions beyond the jel- 
lium sphere lowers the electron density as compared to 
the bulk value 3 -^ (see Fig. [51 inset). On the other hand, 
assuming the tilt (j 1 T[) and not considering density gradi- 
ents in the energy functional may become problematic. 

A simple way to approximately include spill-out effects 
is to use an electronic density which is slightly lower 
than n. In the case of the surface plasmon, where the 
numerically obtained frequency uj sp is lower than u>m, 
such an approach leads to a reduced frequency ujm — 
wm\T — N on t/N, where A out is the number of electrons 
outside the jellium sphere. However, o; sp is still lower 
than ujm , and moreover it exhibits a non- monotonous be- 
havior not accounted for by uim (see Fig. [2] and Ref. \24 ). 
Assuming that the effect of spill-out on the spin-dipole 
frequency is similar to the one on the surface plasmon 
frequency, it is tempting to substitute ujm by u> sp in Eq. 
I2TH As shown in Fig. [21 such an approach (circles) con- 
siderably improves the estimation of w s( j. 
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B. Non-uniform ground-state density 

A further improvement of the accuracy can be achieved 
by going beyond the approximation of the tilt (fT7|) of 
the spin up and down densities and, at the same time, 
taking into account the spatial variations of the ground- 
state electron density. The latter consideration is crucial 
since in the spill-out region the density falls rapidly to 
zero (see inset of Fig. such that the kinetic energy 
contribution Ek.g of (|A2[) . which includes the gradients 
of the electronic densities, becomes important. 

In this section, we assume that the magnetization pro- 
file of the surface paramagnon is given by the static mag- 
netization induced by a static external dipolar magnetic 
field. Expressing the energy functional of Eq. lAll in terms 
of the charge and magnetization densities, the ground- 
state conditions for n (r) and mo(r) are 



5E[n, m] 



Sn(r) 



5E[n, m] 



=n (r) 
=m (r) 



Sm(r) 



=n (r) 
=m (r) 



0. 



(21) 



Applying an external magnetic field along the z axis, 
B ex (r) = £? ox (r)z results (with gn^ = 1 in atomic units 
and the negative charge of the electron) in an additional 
contribution to the total energy functional 



ET[n,m] — E[n,m] + / dr B cx (r) m(r) 



(22) 



If we work with a spherically symmetric nanoparti- 
cle, the charge and spin channels are decoupled. Thus, 
in linear response, the application of a magnetic field 
does not affect no(r), and we drop this functional vari- 
able hereafter. The magnetization density is driven from 



its ground-state value mo(r) = to a perturbed value 
Am(r), which is given by 



6E[r 



5m(r) 



--Bex(r) • 



(23) 



m— Am(r) 



Once the applied field is removed, the nanoparticle is 
left with an extra-energy 



AE 



I [a [a' 

2j r J r 5m(r)6m(r>) 
- J dr B in (r)Am(r) . 



Am(r)Am(r') 



(24) 



In the last equality we have used the perturbed equilib- 
rium condition (|23|) and defined B ln = —B cyi as an inter- 
nal field that counterbalances the applied one. Once the 
perturbation is switched off, the dynamics of the magne- 
tization is determined by the excess energy AE. 

For the dipolar excitations that we are interested in 
{i.e. Eq. ll2b[) an appropriate choice for the external field 
is .B ex = —zjX Bl with 1/Xb measuring the strength of 
the perturbation. With (JT5J) , this allows us to write the 
extra energy as 



AE = — — [ dr z Amir) = 
2Ap J w 



ZAb 

The condition (|23| and the form of _B CX fix the 
induced magnetization, which we can write as 
Am(r) = 6mn(r) cos 9. Defining the quantity m(r) = 
Asr 2 <5mri(r), we show in App.lBlthat in the linear regime 
the magnetization profile is determined by the differential 
equation 



D(r) rh'(r) + 



D(r) ( - + 36D(r) ) + A KS {r) 



where 



A K s(r) 
D(r) 



1 \ 1/3 1 



12tt 



3^\ 2/3 1 



1 n' (r) 
36 n (r) 



m{r) = no(r)r 3 , 



l + c(r.(r)) 



(26) 

(27a) 
(27b) 



with no(r) being the electron density in the ground state 
which can be calculated numerically from a static LSDA 
code. The primes denote derivatives with respect to r. 

As in the simpler case of a tilt-like magnetization, we 
assume that the functional form of the magnetization 
profile is conserved, up to an overall factor, in the os- 
cillations occurring when the external field is switched 
off. Such an assumption is supported by the numerical 



results shown in Fig. [3] The magnetization profile ob- 
tained at resonance (uj — cj s d, filled circles, right scale) 
is much stronger but very close in shape to the static 
one (u> = 0, i?cx = —z/Xb, empty circles, left scale). It 
is important to notice that the solution m(r) of Eq. [^U 
(thin solid line) is a good representation of the local spin- 
density calculations. The differences for small values of 
r are not significant because of the volume integrals that 
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FIG. 3: Radial part of the magnetization profile for vari- 
ous particle sizes. The filled circles represent the TDLSDA 
calculations near the frequency of the surface paramagnon 
(right scale). The empty circles are the static results ob- 
tained from LSDA with an external field _B CX = — z/Xb (left 
scale). The solid line is the solution of (J26J based on a non- 
uniform ground-state electron density. The normalization of 
the magnetization profile differs from that of Fig. [T] by a fac- 
tor (r/a) 3 , such that the tilted spin densities fT7{ result in a 
step function. 



are performed. In addition, we see from Fig. [3] that the 
various approximations for the magnetization profile do 
not deviate considerably from the simple tilt (fT7| that 
we used in the previous chapter. Even if the magneti- 
zation profile attains its maximum value around r = a, 
the spin dipole is not a surface mode (in contrast to the 
surface plasmon), since the excitation is not confined to 
the surface but appears in the whole nanoparticle. 

The restoring force associated with AE, Eq. fJH leads 
to oscillations of the two spin populations with a fre- 
quency 



2 AE 

N ZT 2 



N 



47r J dr rm(r) 



(28) 



Using the profile rh{r) from (|26|) we obtain a good ap- 
proximation (pluses in Fig. [2]) of the numerically obtained 
Lj s d- This shows the importance of the spill-out in deter- 
mining the frequency of the spin-dipole excitations, and 
underlines that the corresponding shift can be accurately 
estimated from the equilibrium density profiles 25 . 



IV. SPIN-DIPOLE SPECTRUM AND 
PARTICLE-HOLE EXCITATIONS 

In the preceding section, we have identified the be- 
havior of the lowest frequency peak in the spin-dipole 
absorption spectrum. Two important questions deserve 
to be addressed now. The first concerns the specificity of 
the lowest energy peak as compared with the other ones. 
The second question, already treated in the literature^, 
partly in the context of electronic excitations in quantum 
dotsii, is whether or not the spin dipole can be consid- 
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FIG. 4: Spin-dipole absorption spectrum (solid line) and 
particle-hole excitation spectrum (dashed line) for two 
nanoparticle sizes. 



ered to be a collective excitation, as it is the case for the 
surface plasmon. 

In Fig. [1] we saw that in addition to its considerably 
larger oscillator strength, the first peak is peculiar from 
the point of view of the induced magnetization which has 
significant contributions of constant sign, whereas there 
are always important contributions of different sign for 
the other peaks. The assumption of a tilt for the magneti- 
zation profile used in our analytical approach of Sec. lIII A1 
is appropriate to describe a mode without nodes and thus 
allows for an estimate of the frequency of the lowest peak. 
In the semi-analytical model of Sec. IIIIB1 the magneti- 
zation is supposed to be generated by an external mag- 
netic field that is linear in z. The radial components of 
the magnetization profiles obtained numerically are also 
positive except for a very small insignificant region close 
to the center in the largest particles (see Fig. [3]) . In order 
to predict the frequencies of higher-energy peaks in the 
spectrum, one would have to assume more complex mag- 
netization profiles. However, those peaks show no obvi- 
ous regularities in their magnetization profile {e.g. the 
number of nodes does not increase monotonically when 
moving to higher frequencies). One can conclude that the 
first peak really stands out as the only one with an es- 
sentially everywhere-positive magnetization, a property 
that allowed us to construct a rather precise theory for 
the frequency associated with this mode. 

The authors of Ref. 8! concluded that the surface para- 
magnon is a collective excitation, based on the observa- 
tion that the resonance exhausts more than 90 % of the 
total spectral weight. A widely accepted criterion is to 
consider an excitation as collective if it results from the 
superposition of a large number of low-energy particle- 
hole excitations. This is certainly the case of the well- 
studied surface plasmo n 26 i 27 , where the residual inter- 
action (understood in this context as going beyond the 
Hartee-Fock approximation) results in a small perturba- 
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FIG. 5: Spin-dipole absorption spectrum for Na34 and in- 
creasing values of a c corresponding to an increasing impor- 
tance of the interactions. 
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FIG. 6: Charge dipole absorption spectrum for Na34 and dif- 
ferent values of a c . For a c — 1, the surface plasmon frequency 
uj sp ~ 0.8 u>m is recovered. 



tion of most of the particle-hole excitations and the ap- 
pearance of the collective excitation in the high-energy 
sector of the spectrum. The considerably lower energy of 
the spin dipole and its size scaling suggest some impor- 
tant differences in the nature of the two excitations. In 
order to test this conjecture, we show in Fig. 2] a com- 
parison between the full spin-dipole absorption spectrum 
(solid lines) and the corresponding particle-hole excita- 
tion spectrum (dashed). The latter is obtained by remov- 
ing the electron-electron interaction in the calculation of 
the linear response (setting a c = in Eq. [7]) , although it 
is still included when computing the ground state. Since 
there is no Hartree contribution for the spin modes, the 
residual interaction in this context is understood as the 
effective exchange-correlation term (JS) . We can see that 
the two spectra have similar structure, with a one-to- 
one correspondence between the excitations. The full 
spin-dipole absorption spectrum appears to be slightly 
red-shifted as compared to the particle-hole spectrum be- 
cause of the attractive nature of the exchange-correlation 
interaction. 

The selection rules for dipole-created electron-hole ex- 
citations dictate a minimal absorption energjsSI, associ- 
ated with a frequency w™ h in ~ (7r/2)(97r/4) 1 / 3 ^- 1 / 3 f s ^ 2 
which has the observed size scaling of cu s d- This estima- 
tion of the first peak of the non-interacting absorption 
spectrum (dashed lines in Fig. [3J agrees within 20 % with 
the frequency that is obtained from Eq. [^D] by only keep- 
ing the one-body (kinetic energy) component. Once we 
consider exchange and correlation corrections, the previ- 
ous excitation splits according to its total spin. The spin 
selection rules tell us that the lower frequency appears 
in the absorption spectrum (solid lines in Fig. [4]). The 
corresponding shift can in principle be extracted within 
the local density functional approximation provided the 
single-particle wave-functions are known. 



In order to investigate in more detail the evolution of 
the single-particle excitations into spin modes we vary a c 
in Eq.[7]from the noninteracting particle- hole case a c — 
to full spin-dipole excitation a c = 1. Fig. [5] shows the 
evolution of the spin-dipole absorption spectrum with a c 
for Na34 . One can see that the structure of the spectrum 
is not modified by the interaction strength. The lowest- 
frequency peak is always the dominating one, and its 
strength is hardly changed. The whole spectrum is red- 
shifted by the interaction, which is globally attractive in 
this case. 

The behavior of the spin dipole can be contrasted with 
that of the charge dipole (Fig. [6]). When a c = 0, the four 
particle-hole excitation modes previously obtained can be 
observed at frequencies lower than 0.5 ~ 1-7 eV. By 
gradually increasing a c , three of the modes are slightly 
blue-shifted and become considerably weaker, while the 
fourth one experiences a much larger blue shift and dom- 
inates the other peaks for a c > 0.2, eventually by several 
orders of magnitude. For a c = 1, this peak coincides 
with the surface plasmon, which is a collective excitation 
with Usp — 2.8 eV « 0.8 wm- 

These findings are at odds with the claims of Serra et 
al&, who interpreted the surface paramagnon as a collec- 
tive excitation. In contrast, our results indicate that the 
various spin-dipole modes appearing in the absorption 
spectrum should be viewed as individual particle-hole ex- 
citations, slightly modified by the electron-electron inter- 
action. 



V. OPEN-SHELL SYSTEMS 

For closed-shell systems, an electric dipole field only 
couples to the charge dipole in the linear regime and the 
excitation of the surface plasmon dominates the absorp- 
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FIG. 7: Dipole absorption spectrum S n (solid line) and dipole- 
induced spin-dipole absorption spectrum Smn (dashed line) 
for Na27. The frequencies of peaks (a), (b), and (c) are 1.22, 
1.82, and 2.77 eV, respectively. 

tion of laser light. One possibility to observe cross-talking 
between the charge and dipole modes is to operate in the 
nonlinear regime with strong excitations^. This may 
raise some practical difficulties, such as electrons escap- 
ing the nanoparticle, thus leaving behind a net positive 
charge. Another possibility is to work with open-shell 
systems that possess an intrinsic magnetization in the 
ground state, so that the charge and spin modes are cou- 
pled even in the linear regim e 12 ' 17 . 

In Fig. [TJ we show the dipole absorption spectrum S n 
and dipole-induced spin-dipole absorption spectrum S mn 
for a Na27 nanoparticle 2 -. In both cases, the system is ex- 
cited by an oscillating electric field which induces both a 
charge-dipole mode (solid curve) and a spin-dipole mode 
(dotted curve). The charge and magnetization profiles 
corresponding to some of the observed peaks are plotted 
in Fig. M 

The strongest coupling between the charge and spin 
channels occurs at the frequency of the surface plasmon. 
For high-energy peaks, the charge-dipole response dom- 
inates its spin-dipole counterpart, while the spin-dipole 
response is more important in the low-energy spectrum. 
The peak labeled (b) in Fig. [7J appears as a special case - 
with an energy intermediate between that of the surface 
paramagnon and surface plasmon modes - for which the 
spin-dipole response is comparable to the charge-dipole 
response. The same qualitative features have been ob- 
served in the spectra of other open-shell systems. 

In the sequel we focus on the maximum of the spin- 
dipole response that occurs at the surface plasmon fre- 
quency for systems with non-zero ground state magne- 
tization (peak (c) in Fig. [7J. For weak excitations, the 
dynamical magnetization can be written as 

Am(r) = A£(r)n (r) + An(r)£ (r) (29) 

in terms of the ground-state electron density no(r), 
the ground-state polarization £o( r )j and their dynami- 
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FIG. 8: Radial part of the dynamical charge (solid lines) and 
magnetization densities (dashed lines) corresponding to the 
peaks highlighted in Fig. [7J The normalizations are as in Fig. 
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FIG. 9: Ratio of the spin-dipole to charge-dipole absorption 
cross-section for a spin-independent excitation at the surface 
plasmon frequency, as a function of the mean ground-state 
spin polarization £. Squares refer to particles Nai3 - Nai7, 
triangles to particles Na2i - Na27, and dots to particles Na4i 
- Na49, where the open electronic shells are d, f, and g, re- 
spectively. 

cal counterparts An(r) and A£(r). As the laser light 
couples essentially to the charge degrees of freedom, the 
excitation has the form of a shift Z of the entire elec- 
tron population. The dynamical excitation density is 
thus concentrated at the surface and, in the hard-wall 
homogeneous density approximation, given by An(r) = 
nZS(a — r) cos 8, just like in the case of the surface plas- 
mon. The resulting charge excitation corresponds to a 
peak at the surface plasmon frequency w sp in the absorp- 
tion spectrum of Fig. [Jj 

In addition, even though the excitation does not act 
directly on the spin polarization (i.e. A£(r) = 0), the dy- 
namical magnetization corresponding to the charge dis- 
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placement Am(r) = An(r)£o(r) does not vanish when 
£o 7^ 0. Assuming that the polarization £o( r ) is uniform 
inside the particle, and equal to the mean spin polariza- 
tion^ = (iV T - N l )/N (for the example of Na 2 ?2£ one 
has £ = 7/27 s» 0.26), we obtain a response to the spin- 
independent excitation I4 x ,n in the magnetization chan- 
nel Am(r) which is proportional to the induced An(r) 
in the charge channel, with proportionality constant £. 
Therefore, we expect 

iaKp) = (S n (u sp ) . (30) 

In Fig. [5]we present results of TDLSDA calculations with 
spin-independent excitation for the mode with lu = lu sp , 
for a variety of open-shell nanoparticles. The ratio be- 
tween the spin-dipole and charge-dipole absorption cross- 
sections for not too large polarization is indeed to a 
very good approximation given by the mean ground-state 
polarization £. This allows us to predict particularly 
strong cross-talk between the spin-dipole and the sur- 
face plasmon modes for open-shell nanoparticles having 
large ground-state polarizations. 

VI. CONCLUSION 

In this work, we have studied the spin-dependent linear 
response in alkali-metal (particularly sodium) nanoparti- 
cles. Our primary aim was to achieve some insight into 
the nature of these modes, which were first investigated 
by Serra and collaborators^. Towards this goal we de- 
rived simple analytical and semi-analytical models that 
were confronted with linear response TDLSDA calcula- 
tions. 

The spin-dipole absorption spectrum displays a num- 
ber of peaks at frequencies lower than the surface plas- 
mon frequency. The lowest of them is characterized by a 
magnetization profile without nodes. An excess of spin- 
up electrons is built in half of the nanoparticle at the 
expense of the spin-down electrons, which are majority 
in the other half of the nanoparticle. The restoring force 
of such a non-equilibrium configuration results in the out- 
of-phase oscillation of the two spin subsystems. The local 
spin-density approximation can be used to estimate the 
restoring force, and within a classical picture, we could 
estimate the lowest frequency. This approach provides 
the correct scaling of the frequency with the particle 
size (as N^ 1 / 3 ), albeit blue-shifted with respect to the 
TDLSDA results. Such a deviation is partially corrected 
by including the spill-out effect in a phenomenological 
way using the numerically obtained surface plasmon fre- 
quencies instead of the Mie value. A more sophisticated 
model, taking into account the inhomogeneities in the 
ground-state density and gradient corrections, yielded an 
even better agreement. 

By comparing the spin-dipole absorption spectrum 
with that obtained by progressively removing the 
electron-electron interaction, we observed a one-to-one 
correspondence of the particle-hole excitations and the 



spin-dipole modes. We thus showed that the spin-dipole 
modes are slight perturbations of the particle-hole exci- 
tations, and therefore do not qualify as genuine collective 
excitations, contrary to the claim of Ref. 8. 

Finally, we studied the possibility of exciting the spin- 
dipole modes by ordinary optical means (laser pulses). 
For open-shell systems, it is well-known that the spin 
and charge modes are coupled in the linear regim e 12 ' 1 . 
We showed that, when exciting the system with a den- 
sity shift, a spin-dipole mode appears at the surface plas- 
mon frequency, together with the standard charge-dipole 
mode. The ratio between the strengths of the absorp- 
tion peak for the spin dipole and charge dipole modes 
was shown to be given by the spin polarization of the 
ground-state. 

While our numerical calculations were done in the case 
of Na nanoparticles, our general conclusions are also valid 
for noble-metal nanoparticles and semiconductor quan- 
tum dots. The latter systems are more adapted than al- 
kaline nanoparticles for experimental spectroscopic stud- 
ies. Moreover, the concepts developed for the study of 
spin modes in normal-metal nanoparticles could be use- 
ful in analyzing ferromagnetic nanoparticles, in view of 
the strong internal field existing in these materials. 
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APPENDIX A: LSDA PARAMETRIZATION 

Our numerical and analytical approaches to obtain the 
spin-dipole resonances are based on the local spin density 
approximation^. For completeness we present in this 
appendix the particular parametrization that we chose 
in our approaches. The energy functional of the electron 
system can be written as 

E[n\n l ] =E K [n r ,n l ]+E K [n}+Exc[n [ ,n 1 }, (Al) 

where Ek represents the kinetic energy, En the Hartree 
contribution, and Exc the exchange-correlation term. 
The kinetic energy is given by 

E K [n\nt] = £ K ,TF[nW] + £K,G[™W], (A2) 
with the Thomas-Fermi component 

Sk.tf^W] = ^ (6tt 2 ) 2/3 J dr (n^ 3 (r) + r^ 5 / 3 (r)) 

(A3) 

and the gradient correction 

^ W] .y»(w + wy (A4) 
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which takes into account the spatial density variations 
in the electron ga o 3 i 29 . This term is particularly relevant 
when we consider spill-out effects and thus a finite region 
where the density exhibits large spatial variations (as in 

sec. hub]). 

The Hartree contribution depends only on the total 
density n = n' + n^, and is given by the electrostatic 
potential energy 



„ r t H 1 /"/",, / n(r)n(r') , , , 

g H [nW] = s //drdr' ]J_\ ' (A5) 



The normalized inter-particle distance r s = (47rn/3) -1 / 3 
and the spin polarization £ = (n* — n^)/n are both 
local properties of the electron system. The above 
parametrization of -Exc was used in our analytical ap- 
proaches, as well as in the numerics, since it has been 
proven to provide a good representation for the electron 
densities that we are interested in. 



This term is irrelevant for spin-dipole excitations in 
spherically symmetric nanoparticles since they do not in- 
volve charge displacements. 

The exchange-correlation term can be expressed in 
terms of the exchange and correlation components as 
Exc = Ex + Eq with 



Ex,cWy]= / drr*(r)e x/c (nW) 



(A6) 



APPENDIX B: MAGNETIZATION PROFILE FOR 
NON-UNIFORM GROUND-STATE DENSITIES 



The function e xc = e x + e c determines the exchange- 
correlation potential through Eq. [5] For the exchange 
part one has 

e x (nW) = -\ (n^ + n^ 3 ) /n. (A7) 



For e r we use^i 



e c \n' ,n 



T n i\ - £ U 



with 



£(r.) 
f(0 



(r.) + [£(r B )-£(r,,)]/(0 (A8) 



-0.1423 



1 + 1.0529^+ 0.3334^ 

-0.0843 

1 + 1.3981^ + 0.2611^ 

(l+Q 4 / 3 + (l-g) 4 / 3 -2 
2 4 / 3 - 2 



(A9) 
(A10) 
(All) 



In this appendix we develop the perturbed equilibrium 
condition (|2"3f for a closed-shell nanoparticle in an ex- 
ternal magnetic field and derive the differential equation 
(|2"6f for the magnetization profile. 

Expressing Eq.[53]in terms of the polarization we have 



1 6E[£\ 



n (r) 6£(r) 



(Bl) 



=A£(r 



Similarly as in (JTo 
tives of E K XF) E x 



we can write the functional deriva- 
and Ec, respectively, as 



1 SE KyTF (3tt 2 ) 2/3 



, . ,,, . n /3 (r) ([l + e(r)] 2/3 -[l-C(r)] 2/3 ) , (B2a 

n Q (r) 5£(r) 4 u V J 



1 SE C 
n (r) 6£(r) 



; -[l-e(r)] V3 ) , (B2b) 
(r s (r))-^(r s (r))]/'(£(r)), (B2c) 



where /' stands for the derivative of the function / defined in (|A11[) . Because of large density variations in the 
spill-out region we consider the gradient correction (|A4p for which we get 

J_«^G_^ S \ f , T { \Vntf |VC(r)| 2 M m _. 

n (r) SS(r) ~ 72n (r) ££(r) [J { n (r) + W 1 - ?(r) ) J ' ( ] 
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The first term of the integrand is independent of £(r), and the second can be treated using partial integration yielding 

1 6E K , G _ 1 / e(r)|V£(r)| 2 l / rc (r)Vg(r) \\ 
no(r) 6Z(r) 36 V [1 - £ 2 (r)] 2 n (r) {!-?(*)))' [ ' 



Since Vno(r) = n' (r)r, we have V(no(r)V£(r)) = 
rio(r)(9^(r)/9r) + no(r)V 2 £(r). Assuming that the po- 
larization is small and a smooth function of r verifying 
V 2 £ <C d£(r)/dr we remain in linear order in £ and write 



1 £E K ,c 
n (r) <5£(r) 



i_ r^w aeoo 

36 no(r) <9r 



(B5) 



with Aks defined in (|27a| . Since Eq. lB6l admits solutions 
of the dipolar form, we write A£(r) = <5£r,(t") cos and 
m(r) — Asr 2 n(r)^R,(r), obtaining Eq.[26]for the magne- 
tization profile. The function rh(r) is more appropriate 
than 5^(r) for numerical calculations dealing with small 
electron densities. 



Gathering the various contributions to the energy, Eq. lBl 
becomes 



1 n' (r) <9A£(r) 
36 n (r) dr 



+ A KS £(r) 



■ cos 6 , (B6) 
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